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We consider a single living semi-flexible filament with persistence length £ v in chemical equilibrium 
with a solution of free monomers at fixed monomer chemical potential pi and fixed temperature T. 

While one end of the filament is chemically active with single monomer (de)polymerization steps, 
the other end is grafted normally to a rigid wall to mimic.k a rigid network from which the filament 
under consideration emerges. A second rigid wall, parallel to the grafting wall, is fixed at distance 
L « £ p from the filament seed. In supercritical conditions where monomer density pi is higher 
than the critical density p\ c , the filament tends to polymerize and impinges onto the second surface 
which, in suitable conditions (non-escaping filament regime) stops the filament growth. We first 
establish the grand-potential fi(pi,T, L) of this system treated as an ideal reactive mixture and 
derive some general properties, in particular the filament size distribution and the force exerted 
by the living filament on the obstacle wall. We apply this formalism to the semi-flexible, living, 
discrete Wormlike chain (d-WLC) model with step size d and persistence length £ p , hitting a hard 
wall. Explicit properties require the computation of the mean force fi(L) exerted by the wall at 
L and associated potential /;(L) = —dWi(L)/dL on a filament of fixed size i. By original Monte- 
Carlo calculations for few filament lengths in a wide range of compression, we justify the use of 
the weak bending universal expressions of Gholami et al.(Phys.Rcv.E. 74,(2006), 041803) over the 
whole non escaping filament regime. For a filament of size i with contour length L c = (i — 1 )d, 
this universal form is rapidly growing from zero (non compression state) to the buckling value 

fb(L c ,£ p ) = ™ ^ over a compression range much narrower than the size d of a monomer. 
Employing this universal form for living filaments, we find that the average force exerted by a living 
filament on a wall at distance L is in practice L independent and very close to the value of the 
stalling force Ff 1 = ( ksT/d ) ln(pi) predicted by Hill, this expression being strictly valid in the rigid 
filament limit. The average filament force results from the product of the cumulative size fraction 
x = x(L, £ p , pi), where the filament is in contact with the wall, times the buckling force on a filament 
of size L c at L, namely Ff 1 = xfb(L-,£ p ). The observed L independence of Ff 1 implies that x oc L~ 2 
for given {£ P ,pi) and x oc In pi for given ( £ P ,L ). At fixed (L,pi), one also has x oc £ p 2 which 
indicates that the rigid filament limit £ p —> oo is a singular limit in which an infinite force has zero 
weight. Finally we derive the physically relevant threshold for filament escaping in the case of actin 
filaments. 
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I. INTRODUCTION. 


Cytoskeleton actin filaments, with the help of a wide variety of auxiliary proteins, are at the root of dynamical 
processes involved in cell motility [H [5]. The growth of lamellipodium and filopedia is directly related to actin filaments 
pushing or sometimes pulling (with the help of trans membrane proteins) with their barbed end pointing against the 
cellular membrane. A subtle interplay of polymerizing or depolymerizing steps, involving single G-actin monomers at 
the barbed end, provides the essential mechanism allowing the cytoskeletal network to keep contact while maintaining 
a permanent pressure force on a load resisting membrane. 

In vitro experiments on biofilaments, like actin and tubulin, to measure in supercritical conditions either the force- 
velocity relationship in detectable non-zero velocity conditions [3] 2], or the approach to stalling where the applied load 
effectively stops the net polymerization of living filaments and the stalling force is effectively measured [5], have been 
deviced. To simplify the analysis and to concentrate on the fundamental process of force generation by polymerizing 
filaments, the experiments deal with bundles of parallel filaments hitting an orthogonal moving wall, a network having 
strong analogy with the structure of actin filopedia [T]. Data analysis requires models for bundle dynamics and stalling 
force predictions, and in general, most models treat living filaments as perfectly rigid. 

In a series of pioneering papers on this topics in the eighties, Hill was the first to propose the expression [B] 

F b H un = N f k ^\np 1 ( 1 ) 

for the force generated by a bundle of Nf growing (proto)filaments stopped by a normal wall. In Eq. ([l]) , T is the 
absolute temperature, d is the effective monomer size along the filament contour (equal to half the G-actin diameter 
as there are two interwined protofilaments in F-actin) and pi = Uq/Wq > 1 is the reduced free monomer density 
equal to the ratio of bulk polymerizing and depolymerizing rates Uq and Wq. Using a combination of thermodynamic 
and mean-field arguments |i.6j, this expression has been established as the equilibrium state (zero growth velocity) of 
a more general expression linking the wall velocity and the load force for a bundle of filaments in a generally non 
equilibrium framework. As stressed by Hill, eq. ([!]) is derived from a one-dimensional longitudinal and incompressible 
model which implies a proportionality between the average polymerization rate and the average wall velocity. 

The ID brownian ratchet model for an individual rigid filament hitting a moving wall was later proposed to offer 
a physically justified stochastic model [7;. The living filament is subject to random polymerizing and depolymerizing 
steps with respective rates Uo and Wo with Uq > Wq to treat supercritical conditions. While the depolymerizing step 
is possible even in presence of the wall, the polymerizing step is only accepted if it does not lead to an overlap with the 
moving wall. In addition the wall undergoes a ID brownian motion characterized by a diffusion coefficient and by a 
load which biases the wall dynamics towards the filament’s end. The coupling between the filament (de)polymerization 
dynamics and the wall random motion leads to a stalling force in agreement with Eq.Q and to a stationary drift 
velocity of the wall, which for large wall diffusion coefficient, agrees with Hill’s prediction of the load-velocity law. 
When many parallel filaments act together as a bundle, the brownian rachet model can be generalized to a multi-rigid 
filament system while remaining essentially ID giiEHini- The dynamical coupling among filaments via the common 
wall evolution, which is very sensitive to the relative longitudinal disposition of the filaments, has strong implication on 
the velocity-load relationship [21 SI El- Let us note that all the above models consider non interacting filaments and a 
single kind of actin-Adenosine triphosphate (actin-ATP) complex for the monomers whether free or incorporated into 
filaments. The above dynamical models can be generalized to take into account lateral interactions among filaments 
mm mainly to treat a many-protofilaments model and/or the hydrolysis of the ATP(Guanosine triphosphate-GTP) 
in the filament actin-complexes (tubulin-complexes) by considering additional types of complexes, requiring in turn 
additional information on specific rates mm- 

If rigid filament models are certainly satisfactory as long as the elementary working filaments (being isolated and 
uncrosslinked) remain sufficiently short, the flexibility of F-actin should be properly considered for longer filaments. 
Flexibility was found to be relevant in some important experimental cases. The bending shape of single F-actin 
filaments observed by fluorescence spectroscopy, was precisely exploited to measure for the first time the typical 
polymerization force generated by single living actin filaments na ns]- In the optical trap experiment of Footer et 
al.[5]j a bundle of about ten filaments, with their seeds glued to a trapped colloidal particle, push with their active 
side (barbed end) on a fixed rigid wall. The polymerization force they progressively develop to reach equilibrium is 
inferred by measuring the colloidal displacement in the trap. The interpretation of the experiments was possible only 
by assuming the presence of escaping filaments in the bundle (filaments growing parallel to the obstacle wall after a 
large angle bending fluctuation), a phenomenon interpreted by the authors as rod buckling related the beam elastic 
instability [I]. Despite care in eliminating data potentially polluted by escaping filaments, the measured stalling force 
for a eight filament bundle was (repeatedly) found to be close to Eq. ([T]) with Nf ~ 1 instead of the expected Nf = 8 
filaments number, a result still presently not understood. 
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That flexibility leads, in some extreme cases, to escaping filaments was reported and analyzed in a non equilibrium 
simulation of a model of single living filament hitting a moving wall in which filament flexibility was explicitly taken 
into account [T7J HB]. Quite generally, in these pseudo-stationary simulations with constant load, a wall velocity 
enhancement was found with respect to the predictions of the "rigid filament-hard wall" ratchet model, in agreement 
with theoretical considerations which have generally predicted an enhancement of the efficiency of the conversion of 
chemical free energy into useful work when realistic filament flexibility is included |1 91120j . For large loads (still below 
the stalling force) and for large seed-wall distances, some escaping filaments were detected during the drift of the 
wall na esi- It was argued that this phenomenon is related but distinct from rod buckling and hence was denoted 
as the "pushing catastrophe". The consensus seems to be that to efficiently grow against membrane resistance, actin 
filaments should be neither too short (short filaments are too rigid to intercalate easily a polymerizing monomer 
between the tip of the filament and the wall) nor too long as the load would simply buckle them, the optimal range 
70 nm — 500?rm being cited in a recent review article [2T:. 

In this paper we concentrate on the equilibrium Statistical Mechanical treatment of a semi-flexible filament in a 
slab. In section [IIJ extending previous work [22], we establish within the reactive grand canonical ensemble, the 
grand potential for a living filament in contact with an obstacle wall at fixed temperature and fixed free monomer 
chemical potential. Formal expressions for the size distribution and the equilibrium force on the obstacle wall are 
established. Section III deals specifically with F-actin modeled as a living discrete Wormlike Chain (d-WLC). We 
first define the model and the related range of physical parameters to probe the non-escaping regime of the filament. 
We then compute, by Monte Carlo simulation, the compression-force law for a dead (non-reacting) d-WLC in the 
slab and validate, in the non-escaping regime, the weak bending expression of Gholami et al [25j . Subsequently we 
define the filament force averaged over a distance equal to a monomer size d , crucial for the comparison with Hill’s 
prediction, In section IV we introduce the stalling force and compare the predictions for flexible (finite t p ) against 
rigid (£ p — > oo) models, proving for the latter Hill’s expression for the stalling force. In the entire range of filament 
lengths corresponding to the non-escaping regime, the flexible filament has a stalling force only few percents larger 
than a rigid filament (Hill’s law). Nonetheless, the specific L-dependence of the force (~ L -2 ) resulting from buckled 
filaments hitting the obstacle wall, induces a spectacular, previously undetected, effect of flexibility. Since the stalling 
force is nearly independent of the slab’s width L in the non-escaping regime, this requires a systematic ev olut ion with 
L 2 of the fraction of sizes of the filament touching the wall. This is discussed in section [VJ Finally, section |vT| provides 
some general conclusions and perspectives on the flexibility issue for many filaments bundles, including both static 
properties and dynamic aspects linked directly to the exploration of the force-load relationship. 


II. THE SINGLE GRAFTED LIVING FILAMENT IN A SLAB SYSTEM 
A. The single grafted living filament concept 

We consider a reacting ideal mixture in a slab at temperature T consisting of N t monomers which can either be free 
(G actin-ATP complex) or integrated within a single self-assembled filament (F-actin) with fixed persistence length 
t v . In the F-actin case, t v = 5370d and d = 2.7 nm is the effective monomer size in the filament. The filament, with a 
variable size i and associated contour length L c = (i — l)d, is grafted normally to one of the walls of the slab considered 
as an orthorhombic volume of transverse area A and width (wall to wall distance) L « t p . The filament undergoes 
single monomer (de) polymerization events with a polymerization rate Uq = k on pi , proportional to the free monomers 
density pi , and a depolymerization rate Wq = k 0 ff , independent on the free monomer density, where k 0n / 0 ff are 
the kinetic constants for the (de)polymerization reactions. Supercritical conditions are realized whenever the bulk 
polymerization rate is larger than the depolymerization rate, which happens for pi > k 0 ff/k 0n = 1/K 0 , where K 0 is 
the bulk reaction equilibrium constant[23]. We define p± = P\Kq as the reduced free monomers density; supercritical 
conditions correspond to p\ > 1. In super-critical conditions where polymerization dominates, the filament will grow 
and hit the opposite wall as soon as L c > L. 

The series of possible chemical reaction will be denoted as 

Ai-i + Ai (2 <i<z*) ( 2 ) 

where Ai and A\ represent respectively the grafted filament of size i and a free monomer. At global equilibrium, 
the chemical potentials pi of the different species involved in any reaction must satisfy the chemical equilibrium 
requirement 

Pi = p^ + pi ( 3 ) 

This series of reactions is considered as limited to a size window going from a minimum filament size of two (to be 
considered as an effective permanent seed of the filament) up to a maximum size z*. Fixing a maximum filament 
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FIG. 1. Two configurations of the same living filament grafted normally to the left wall through its first two monomers drawn 
in dark grey, in chemical equilibrium with a free monomers solution. According to the whole set of (de) polymerizing reactions 
Eqs. ([ 2 ]), the filament can polymerize or depolymerize at its free end by addition or removal of one single monomer as illustrated 
by the arrows. In our illustrations (a,b), the contour length L c of the polymerizing filament in (a) appears to be longer than 
the distance L between the walls while in (b), the contour length L c of the depolymerizing filament is shorter than the distance 
L. 


length is only necessary when considering flexible filaments. In fact for rigid ID filaments, the obstacle hard wall will 
necessarily limit the filament growth. Instead for flexible filaments, equilibrium Statistical Mechanics based on the 
concept of a steady equilibrium state, can only be applied if we consider a mechanism limiting the filament growth in 
particular if supercritical conditions are considered. Again if the slab is narrow enough with respect to the filament 
persistence length, the obstacle wall will effectively limit the filament growth but for wider slabs we need to introduce 
an artificial limit. In our present geometry (see fig. [T]) we impose a maximum filament size 

C = I NT (g) . (4) 

where INT{x ) means the integer part of the argument. In absence of this size limit, corresponding to filaments 
adopting a planar configuration which covers a quarter of a cycle of radius L : filaments longer than z* could grow 

unhindered parallel to obstacle wall, preventing the establishment of an equilibrium state. The imposed upper limit 

will bias the statistical mechanics averages except for suitable external conditions (choice of control variables N t /AL 
and L in particular) for which the statistical weight of filaments of size z* or longer is negligibly small. This defines 
what we call the "non-escaping" regime at stalling conditions. 

In terms of the temperature T, volume V = AL, total number of monomers N t and total number of grafted filaments 
Nf, the reversible change of the relevant Helmholtz thermodynamic potential F R for this reactive system (hence the 
R superscript) is 

dF R = —SdT — pjyAdL + p?LdA + p±dN t + — 2/ii)dNf (5) 

where S is the system entropy and where p^A and ptA are the total normal and tangential forces exerted by the 

obstacle wall on the system. The last two terms involve the chemical potential pi of free monomers and the chemical 
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potential /i 2 of grafted filaments (seeds) of minimum size 2. These two last terms result from imposing chemical 
equilibrium Eqs.(j3]) for all reactions ([ 2 ]) to the original series of terms Xu=i (J-idN involving all species of the mixture. 
In the latter sum, the chemical potential of any filament with size * > 2 is substituted by = H 2 + (i — 2)/ ii which 
follows recursively from Eqs.([3]) and a final regrouping of terms yields Eq.([5]), taking into account the expressions of 
the total number of monomers N t = Ni + X)~ =2 i Ni and of the total number of filaments Nf = Y^i =2 Ni. 


B. Free energy of a grafted living filament under confinement 


Applying equilibrium statistical mechanics to a closed reacting ideal system [23], the canonical partition function 
Q r = exp [— /3F r \ for a single grafted filament in a solution of free monomers is given by 


Q n (A,L,T,N t ,N f = 1) = 


9i 


(N t — 2) 


tQ 2 


<T- 3) 


93 


Qi 


C Nt-z *) 




( 6 ) 


(N t -2)\^ ' (N t - 3)!- 10 ' "(N t - z*)\ 

The sum over all distinct microscopic states compatible with the macroscopic variables is expressed in Eq. © as a 
sum over (z* — 1 ) similar terms, each of them corresponding to one particular size of the single grafted filament and 
the remaining free monomers. Each term of this ideal system involves the canonical partition function 99 (A, T) of the 
filament of size 2 < i < z* grafted in the slab and the corresponding contribution from the free monomers 


qi (L,A,T) = ^ 


(7) 


where A is the free monomer thermal de Broglie wavelength. 

To each term qi(L,T) corresponds a canonical partition functions qf(T) relative to the same grafted filament of size 
i in the absence of the opposite wall. Keeping the temperature dependence implicit, we define the ratio’s 


/• r \ Qi( L ) 

oc(i,L) = — Q- 
Ti 


( 8 ) 


As long as the intra-filament interactions have a local and homogeneous character, t he ra tio between successive parti¬ 
tion functions q®_ 1 and q { - is independent of i. Hence, as further detailed in section 
dependent equilibrium constant Kq [23] 


IID 


9? 


n° 

jy- — \ 3 

^ V ° — n° n /T/ _ n ° ^ 

q.i-iqi/V 9,-1 


we introduce a temperature 

(9) 


Using Eqs.(8|9l, the partition functions of the filaments of any size i can be written as 


Kt 


qi(L) = a(i, L)ql ( ) (2 < i < z*). 


i-2 


( 10 ) 


where 9 ° is the partition function of the grafted seed. 

Eq. ([g]) can now be combined with expressions TO , to give 


Q H (A,L,T,N t ,N f = 1) = q 0 2 q^~ 2 x 


«(2 ,L) , q^a(3,L)(^) q~^ ^ 2) a(z* , A) (^f* 2) 


= q^qi* 2 x 


9°9f* 2 


(Nt — 2)! (JV t -3)! " (N t — z*)\ 

a(2,L) , a(3,L)(fy) a(z*,L) (K?) {z *~ 2) 


(N t ~ 2 )! (N t - 3)! 


(JVt-2)! 


cr(2, A) T o(3, A) 


Ko 

V 


(N t - z*)\ 


(N t — 2) + ..a(z*,L) 


Ko 

V 


(z*-2) 


(N t - 2)! 
(N t - z*y. 


and the partition function can further be transformed as 


0y2 N t -2 

Q*= q l V n , Q } _ x 


Kq (N t -2)1 
q° 2 V 2 q^~ 2 ^ 


Kn 


K 2 N t \ 


a(2,L) ( § 


K, 


( 2 , A) -A +a(3,A) (N t — 2) + ..a(z*,L) —f 


N t \ 


V ) (N t - 2 )! 


V 


+ a(3, A) [ —"y 


KoV (N t -2)1 


V J (N t -z*)\ 


N,\ 


V J (N t - 3)! 


+ ..a(z*,L) ( ^ 


N,\ 


V J (N t - 2 *)! 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 
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In the thermodynamic limit (T.L.), here N t —> oo, A —> oo with fixed ratio N t /A = p t L , one gets 


Q R (A,L,T,N t ,N f = 1) 


g°A 6 
Kq N t \ 


x D{p t ,L) 


(16) 


where we have defined 

D{pt,L) = a(2, L)pf + a(3, L)p ^+ a(z + 1, L)p[~ + 1 - ) + .. + a(z*,L)p[ z •* (17) 

Pt = PtKo = (18) 

and where we have replaced ~ 1 in all terms of D(p t , L). 

We have thus 


f}F K (A, L, T, N t , N f = 1) = N t [In (A 3 p t ) - 1] - In - In D(p t ,L) 


(19) 


Note that in the T.L. we may replace N t by N\ = N t — Ifu (the average length of the filament) and pt by pi, so 
that the first term in the r.h.s. of Eq.( 19 ) is the Helmholtz free energy of the bath of free monomers. While the third 
term of Eq.(19) is the relevant free energy of the living filament, the middle term, function of T only, must be linked 
to the free energy required to graft the filament seed (fixed dimer). 

The probability for the living filament to have a size j, defined as P(j) = P(j ; L , p t ), is the term of index j in the 
global partition function Eq. (J6]) , properly normalized. Using the equivalent version of Eq. (16), leads to 


Pii) = AifA 


( 20 ) 


for j £ [2, 2 *], where D is given by Eq.(|17|. 

Performing a Legendre transform of the reactive Helmholtz free energy F R to the reactive grand potential tt R = 
F R -p x N u Eq.§ becomes 


dCl R = —SdT — p N AdL + p T LdA — N t dpi + (p 2 — 2pl)dNf 


( 21 ) 


To obtain fl R , one needs to express pt\ in terms of the old variables according to Eq. ([ 5 ]) using Eq. (jl9]) for F R and 
associated D and p t Eqs. ( fl7ll8t . We get 


Pt) Ajmm Ko/v 


dN t 
= In (A 3 p t ) - 


D(Pt) 


Dipt) 

-h,(AV0-£ 


PtdD(p t )/dpt 1 ^ A 3 „ \ l fa 


where we have introduced the average length of the filament 

EL 2 i «(J» L ) Pt 


lfilial Pt) 


D 


Formally the Legendre transform requires the inversion of Eq.(23 > as N f = N t (p 1 ) to estimate 
/3 U r (A, L,T,p 1: N f = 1) = [pF R (A,L,T,N t ,N f = l)-»4 (=%i) 

Using Eqs.( 19|23 l, one gets successively 


( 22 ) 

(23) 

(24) 

(25) 


(m R (A,L,T,p u N f = l) 


-N t + lfa - In! ^j - In D{p t ) 


N t =Nt(m) 


-Nt + lfa-- 


N t 




J Nt=N t (Mi) 


(26) 
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where we have developed D(p t ) around p\ up to first order. We also note that Eq. (23) can be rewritten as 

2 


Pp l = In (A 3 pi) + O ~ ln (A 3 pi/A' 0 ) 


(27) 


where the central expression is the expected relationship for the chemical potential of one specie s in an ideal mixture, 

the negligible correction coming from the approximations made earlier to simplify the Q R Eq. (15). 

;2 ri — 

Neglecting the term 4A in Eq. (26), the grand potential can finally be reformulated as 


N t 

iR 


pQ*(A,L,T,p 1 ,N f = 1) = 




V K% 


ln D(pi) 


= exp (Pm) 


(28) 


In the biophysics literature it is customary to use the reduced free monomer density p\ as the independent variable 
instead of the more appropriate chemical potential p,\. Therefore we will re-express the grand potential in Eq.(28l as 

pn R (A, L, T, /ir, Nf = 1) = /30 frae (A, L, T , pi) + pn m (L, T, pi) 


A T 

pQ bee (A,L,T,p 1 ) = -—p 1 
A o 


n a {L,T,h) = -k B T 


In 


fq°2 A e 
{ K 2 o 


ln D(pi) 


(29) 

(30) 

(31) 


where we indentify the grand canonical contributions fl h ' ee (A, L, T, pi) and fl m (L, T, pi) for the free monomers and 
the grafted living filament respectively. 

The filament size distribution, the normalization factor D and the average size of the filament, given respectively 
by Eqs.( 20|17|24 1, take the final form 


P(j)=P(j-L,p i) = 

Z* 

j=2 


Nj _ a(j,L)p{ 


N f 


D 


lfil(L,p l) — J 


EL 2 J a ih L ) p\ * d\nD 51n D 


D 


= Pi 


(32) 

(33) 

(34) 


dpi dPp i 

where p\ in the r.h.s. is again used instead of p\ and where Nj is the average number of filaments with size j within 
the microscopic states of the reactive grand canonical ensemble. 


C. Single filament force exerted on the opposite wall 


Combining Eqs. (211 and (29), and noting that pnA is the sum 
filament average force f±(L, piJTone gets 


of a free monomer contribution and the single 


Pf±(L,p 1 ) 


(d(0Qf a )\ fd\n D\ 


(35) 

V dL ) { dL ) 

D 

Z. dL p (r,L,Pi) = 

Z* 

Y J mL)p(i-Lr P i) 

(36) 


i =2 i—2 


where we have used Eq. (321 and introduced a filament mean force potential and associated mean force at fixed length 


Wi(L) 


ML) 


—fcsTln a(i, L). 
dWi(L) 
dL 


(37) 

(38) 


Eq. ( |36| ) gives the equilibrium force exerted on a living grafted filament by a fixed planar wall located at a distance L 
from the grafting wall. As expected, it is the average of the force exerted by the wall on a fixed length "dead" grafted 
filament (this latter force is an average over its internal degrees of freedom), weighted by the absolute probability 
P(i; L,pi) of having a filament of length i. Of course, only the filaments sufficiently long to interact with the wall 
(a(i,L) ^ 1) contribute to the average. 


























D. Equilibrium constants and rates 


Here, we discuss a few properties of the equilibrium constants valid for the arbitrary grafted flexible filament. 
For the considered ideal mixture of Ni free monomers and the series of Ni grafted filaments of size i, one has [22 


13hi = - In 


<h_ 


= In A 3 


PlH = - In TF =-ln 


f In pi 
a{i, L)q .? 


N- 


f 


In P(i) 


(39) 


(40) 


using Eqs. d7l8pt - 

Substituting Eqs.(39|40) in Eq.([3| for any reaction given by Eq.([2]), one gets 


Ki(L,T) = 


P(i) 


a(i, L) q° A 3 


a(i, L ) 


P{i~P)Pi ol{% 1, L) q^_i a(i-l,L) 


Ko(T) 


(41) 


which defines the equilibrium constants K % and its link with the equilibrium constant /to, already defined in Eq. & 
which would apply in absence of wall. Eq.(41l expresses the evolution of the equilibrium constant Ki(L,T) with 
increasing i, as a result of interferences between filaments of sizes i and (i — 1) and the wall. 

Some considerations on the related wall influence on the (de)polymerisation reaction rates are provided in appendix 
A, given their close connection to the equilibrium constants /Q. These rates become essential ingredients of the present 
approach when extended to the study of the coupling of a mobile wall dynamics and the filament (de) polymerization 
steps. 


III. THE WORMLIKE CHAIN MODEL AND THE F-ACTIN CASE. 

A. The discrete model. 

To model the living grafted filament with fluctuating size in the range 2 < i < z*, we adopt the d-WLC model with 
discrete contour length step d and persistence length i p . Using a cartesian reference frame where the grafting wall is 
at x = 0 and the obstacle wall at x = L, the filament normally grafted at the wall at x = 0 has its two first monomers 
located at f-\ = (0,0,0), = (d, 0,0). The filament with i monomers, having a contour length L c i = (i — 1 )d has a 

configuration fully specified by the set of coordinates \fj\j=i,i including the grafted dimer. The instantaneous internal 
potential energy of the filament of size i is 

£(Hi) = -(* - IK + * £[i - cos6> fc ] (42) 

a k=2 

where k = ksT(. p is the bending modulus of the filament, eg the bonding energy associated to the chemical step 
Eq. ([2]) [221 an d 0k the angle between successive bonds implying monomers (k — 1, k, k + 1). The configuration having 
the minimum energy corresponds to the straight filament with all bending angles at zero. The monomer-wall potential 
is zero or infinite depending whether the articulation point (monomer) j is in the slab space (0 < Xj < L) or lies 
inside the obstacle wall (xj > L). If we represent by U w the global filament-wall interaction potential, being the sum 
of all monomer-wall potentials, according to Eq.([8]) the factors a(i,L) become 

i 

a(i,L) =< exp —[PU W ] > ij0 =< JJ©(L-*_,-) > i)0 (43) 

i=2 

where < ... >i_o denotes a canonical average with weight exp (— f3E{[r\i)) (Eq.(|42|)) of a grafted filament of size i in 
absence of the obstacle wall, and 0 is the Heaviside function. For this model, of contour length L C i = (i — 1 )d, we 
have qi(L) = qm and hence a(i, L) = 1, as long as i < z, where z is the integer given by 

2 = INT(L/d) + 1 (44) 

In the case of a living filament undergoing (de) polymerizing reactions and for short enough filaments (?’ < z), this 
WLC model leads to the following expression for the equilibrium constant, as defined in Eq. & m 

d 4 d 4 

K 0 = 2n exp (/3e' 0 ) — [1 - exp (-2 £ p /d)\ « 2?rexp (^e' 0 ) —. 

•tp Zp 


( 45 ) 
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in terms of the fundamental parameters c?, i p , eg of the filament model. The equilibrium constant for filaments hitting 
the wall is given by Eq. (41) where the a factors are given by Eq.(43). 


B. Explicit calculations for the F-actin case. 

1. The relevant L and pi regime to probe single F-actin polymerization force. 

The essential ingredients to get static properties of grafted actin filaments are the wall factors a(i,L) for a grafted 
d-WLC hitting a hard wall (z* > i > z), with i v = 5370c? and d = 2.7 nm. On this basis, all properties can be derived 
for any supercritical value of the reduced monomer density p\ > 1. 

We first consider the relevant range of wall position L and the range of reduced free monomer concentrations p\ for 
which the polymerization force is operative and of interest for a quantitative comparison with in vitro experiments [1] 
0 i n>]. According to Mogilner [2T], to produce a working force, individual filaments should be longer than L c = 70 nm 
(about 25 monomers) to avoid being too rigid but should remain below L c = 500 nm (about 185 monomers) to avoid 
what he refers to as buckling. In Footer et al. experiments [5], the polymerization force was measured for non buckled 
filaments of length 200nm (about 70 monomers). In filopedia bundles [2T], parallel filaments are cross linked by fascin 
but free portions of filaments at the leading edge are supposed to be of the order of 20 — 200nm. Finally in the recent 
experiment of Demoulin et al.0 the bundle length studied to get the velocity load relationship is of the order of 
100 — 400nm (about 40 — 150 monomers per filament) (see supplemental information of ref. 0). Further, it has to 
be noted that in vitro experiments probe the polymerization force in moderate supercritical conditions (1 < pi < 3) 
to avoid too rapid buckling and interferences with spontaneous nucleation of new filaments 00 - 

In our illustrative section of the F-actin case, we will concentrate on the supercriticality regime by considering two 
values of the reduced density pi = 1.7 and p\ = 2.5. We will be interested to the wall position regime 20c? < L < 100c? 
where actin filaments are sufficiently long to avoid unphysical influence of minimum size filaments (j = 2) but still 
sufficiently short to avoid escaping filaments, as it will be made more precise later. 


2. The compression law of a grafted (fixed size) filament. 


The basic input of our theory are the a functions (see Eq.([8])) of "dead" filaments of contour length L c = (n — l)c? 
(z < n < z*). The force /(L,T; L c ) exerted by a wall located at L on a d-WLC filament with £ p = 5370c? and of 
contour length L c > L (see Eq. ( |38| ) ) (using here a notation without size index as we now deal with a unique dead 
filament hitting the wall) has been computed by Monte-Carlo simulation. The resulting force-compression laws for 
three filament sizes (n = 41,77,158) are shown in Figure [2] The MC sampling was realized by a mixture of two 
types of attempted moves, i) local crankshaft moves, where a sequence of three, four or five articulation points are 
rotated as a rigid body around an axis joining the two surrounding articulation points, and ii) pivot moves implying a 
global rigid rotation around a bond of the end chain fragment starting from that bond (the size of the fragment being 
sampled between 1 and (n — 3) articulation points). The force exerted by the filament on the wall was estimated as 


f{L,T; L c )/k B T = lim - In 
A—>0 A 


q(L,T- : L c ) 
q{L — A, T; L c ) 


(46) 


where g(L,T;L c ) is the partition function of a single grafted filament of contour length L c . This force is easily 
estimated during the MC sampling by measuring the probability that the filament configuration has an articulation 
point located in the region of thickness A adjacent to the wall. 

In Figure [2 we observe that as L decreases down from L = L c (where the force vanishes) the force quickly increases 
up to a pseudo-plateau before undergoing a final steep rise as L approaches a value of « 10c? — 15c? on its way down 
to zero. As we are interested to filament lengths limited to z* = I NT (| L/d) corresponding to the non escaping 
regime, it implies that only the elasticity law for the regime L c < L < ^ L c needs to be exploited for each L c . The 
relevant regime for our present study terminates within the pseudo plateau at L c — L = (1 — 2/7 t)L c indicated by a 
vertical arrow for each of the lengths reported in the Figure. As it will be discussed elsewhere [25], the fast increase 
of the elasticity force at short L corresponds to a L c independent behavior / = AksTl p /L 2 , with A ss 2/3, valid for 
escaped filament lengths L c /d » z*. 

Gholami et al. m derived a weak bending approximation for the elasticity law of a grafted continuous wormlike 
chain. Their prediction for the L « i p regime, leads to the universal law summarized here. Using notations of 
reference (25], the identification with our formalism for a dead filament of size n and contour length L c = (n — 1 )d 
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FIG. 2. Compressional force f(L,T;L c ) in units ksT/d, exerted by a grafted dead d-WLC filament of contour length L c 
(L c = (n — 1 )d with n = 41,77,158) and persistence length l v = 5370d on an obstacle hard wall oriented normally to the 
filament grafting direction and located at distance L from the filament’s seed. For each filament size considered, we observe 
three successive regimes as L decreases from L c . First, a rapid rise corresponding to the weak bending regime, followed by 
a pseudo-plateau regime which terminates with the onset of the escaping filament regime characterized by a 1/L 2 divergent 
behavior as L approaches 0 (as (L c — L) approaches L c on the Figure). The L threshold value where the filament enters the 
escaping filament regime is indicated by a vertical arrow for each filament length (see text). In inset, all force data f(L,T; L c ) 
are reduced by the theoretical weak bending plateau value fb(L, L c ) for the continuous WLC model (see Eq. ( |52[) ) and plotted 
versus the renormalized compression distance r/ (see Eq.( |48| )) in order to test the weak bending expression Eq. ( ]53| ) shown by a 
continuous green line. 


leads to 


a(n,L ) = Z||(t 7 ) 


(47) 


in terms of a new reduced compression distance 


L c -L 



involving the characteristic length 



The central quantity Zy is (see eqs. (36) or (38) in reference [25] ) 


Z\\(v) = 2 [(-!) fc+1 A fe 1 exp [-A 2 k fj}\ 

k=1 


(48) 


(49) 


(50) 


where A* = (2fc — 1)|. The (microscopically) averaged force f(L,T ; L c ), defined by eq.(38), which is the force exerted 
by the wall on a (dead) grafted WLC filament of contour length L c hitting the normal hard wall at seed-wall distance 
L, is[25] 


f(L ; L c ) = / 6 /|| (fj) 


( 51 ) 
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where fb turns out to be equivalent to the buckling force for a clamped rod of contour length L c jTj, namely 


fb = 


it 2 ksT 


(52) 


and where f» (fj) is a universal function defined by 


~ 4 d\n.Z\\{rj) 


(53) 


This function, shown in the inset of Figure [2] starts from 0 at fj = 0 and increases monotonically to a unity plateau 
which is reached around fj = 0.25. We argue that the Gholami et al. elasticity function is worth exploiting not only 
for the weak bending regime (limited to fj ~ 0 . 6 ) where it is rather precise, but also for the intermediate pseudo 
plateau regime up to L c /L = 7 r /2 ss z* / z where the force appears to be underestimated by 10 — 15 percent only. 
When this approximation is made for our purpose, the gain is enormous as we do not have to run a large number of 
single filament MC simulations to get the force for each specific filament size n ( L c = (n — l)d) as a function of the 
continuous L variable. We get all the needed expressions as functions of a single universal variable fj (Eq.(48l) under 
the form of an explicit convergent series easy to compute (see Eqs.(50|53 1). 


3. Living filament force in the grand canonical ensemble and the L average force concept. 


Given the properties of the WLC discussed in the the previous subsection, the general expression of the polymer¬ 
ization force for the d-WLC model, Eq. (36), takes the explicit form 


.2* 

/3/j_(L,pi) = Pf^ L ) P & L iP i) 

i=z+l 


(54) 


where z is defined in Eq.(44l and where ji(L) and P(i~,L,pi) are given by Eqs.( 38|32|33 1, computed with model 
functions Eqs.( 51|47 l. 

In Figure [3] w e report the polymerization force f± for the d-WLC model adapted to F-actin persistence length, as 
given by Eq.( |54| ), for two values of pi and highlighted within specific ranges of L. We observe a pseudo periodic signal 
of period A ~ d. In the lower L regime (upper panel), only the first term in the rhs of Eq.( |54| ) (i = z + 1) contributes 
to the force. To discuss the behavior of the polymerization force, let us focus on the interval where L changes from 
L = 28 d up to L = 29d (where 2 = 29), noting that the force f±(L) is essentially zero at both boundaries. In this 
interval, f±(L ) oc f^(L)a(2>9\ L)/D(L). Given the moderate variation of D with L and the fact that the force / 30 
is essentially constant except when L approaches the filament contour length L c = 29d by less than 0.05d (fj ~ 0.3) 
where it drops quickly to zero when L reaches L c , the rise of f±(L) reflects the L dependence of P(30; L, p\) and thus 
the dependence of a(30 ,L) from practically zero (for L < 28.5 d) towards unity. The fast drop of f±(L) towards zero 
at L = 29 d comes from the drop of / 30 . As the L domain increases (lower panel), the variation becomes more complex 
as it involves increasingly more terms in Eq. (54). We note in Figure [3] that Hill’s prediction for the polymerisation 
force lies very close to the average value f±(L) within any interval. 

In his seminal paper, Hill[ 6 ] introduced the stalling force through the work needed to add reversibly a new monomer 
to a rigid filament pressing normally against a wall, as the wall moves by a distance d (See Eq.2 in refjB]). This 
implicitly defines the average of the living filament force over an interval [L, L + d]. On the basis of a reversible change 
of the grand potential Eq.(21 1 , the reversible work at constant T and constant pi performed by the filament pressing 
against a wall moving from L up to L + d satisfies 


W L , L+d (L , pi,T) = n pl (L, pi,T) - n pl (L + d,pi,T) = k B T In 


D(L + d , pi) 


D(L,pi) 


= F%d. 


(55) 


where we have defined the average force over the d interval. 

The average force concept appearing in Eq. (55) can be used for rigid filaments but also for flexible filaments as 


modeled by the d-WLC model, if one adapts specifically the D(L, pi) terms (in particular the a functions) in Eq.(33). 
In order to do so, we first reexpress this average force in an equivalent but more appropriate way to discuss the 
specificities of filament flexibility. Given the definition of z in Eq.(44l, we have z(L + d) = z(L) + 1. For the upper 
limit z* given by Eq.Q, one has z*(L + d) = z*(L) + k where k = 1 or k = 2 depending upon the precise L value 
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L/d 




FIG. 3. L dependence of the polymerization force Eq.(36l for a F-actin single living filament, modeled by the living d-WLC. 
The force is shown for two values of the reduced free monomer concentration, pi = 1.7 (black curve) and pi = 2.5 (red curve), 
as a function of L in the interval [27d, 30d] (upper panel) and in the interval [67d, 70d] (lower panel). The dashed horizontal 
lines show Hill’s predicted value for the reduced force In pi at each considered free monomer reduced concentration. 


which is truncated by the integer value operator in eq. Q. We can then rewrite D(L + d ) in alternative equivalent 
ways: 


(56) 

(57) 

3 (**(L+d)-1) ( 58 ) 


Dropping the last term, present only when z*(L + d) — z*{L) = 2 but which is anyway negligibly small in non escaping 
regime conditions, it leads to 



z(L-\-d) 

z* (L-\-d) 

D(L + d) — pi 

E « _1 + 

^2 a(j,L + d) p-p 1 


*= 2 j- 

—z(L-\-d)-\-l 


z(L) 

z* (L-\-d) — l 

= Pi 

Pi + E p? + 

a(n + l,L + d)p 


m—2 

n=z(L )-\-1 


z(L) 

2*(L) 

= Pi 

pi + Y,pt + 

y a(n+l,L + d)pi 


m—2 

n=z(L )-\-1 


D(L + d) — pi 


z{L) z”{L) 

Pi + ^ PT + a(n+ 1,L + d)p" 


m—2 


%=z(L)-\-1 


(59) 


PdFft{L, pi,T) = In p\ + In [pl + 


( 60 ) 
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where use of Eq.(55 > has been made and where D shl P is obtained by substitution of all terms a(j,L) in D by 
a s hift (j,L) = a(j +1 ,L + d) for all j > z(L). Note that the first term in eq. ( |60| is the single filament polymerization 
force of Hill. 

A similar calculation for the variation of the average size of the filament as the wall moves reversibly from L to 
L + d gives, according to Eqs.(34|55|60), 


pi) — lfu{L d, pi) lfu(L, pi) — 


d 




■ In 


D{L + d, pi) 


D(L,pi) 


d 




-[f3WL,L+d\ 


= 1 + Pi wz- In 

dpi 


\pi + D shift ( L ,Pi)} 

D(L,pi) 


(61) 

(62) 


where the first unity term is also the Hill’s result for rigid filaments hitting normally the obstacle wall[6]. 

In fact, in both Eqs.(60l and (62), the second term on the r.h.s. gives the correction arising from two different 


effects. The first effect is the imposed minimal size of filaments which manifests itself at low L by the additive term 
pi in the numerator of the argument of the logarithm. The second effect linked to the ratio — d(L pi) the effect 

of flexibility by opposition to the purely rigid case (£ p = oo) where D shl P = D. Before embarking on this analysis, 
presented in Section [IV] we derive the precise criteria to be satisfied, in supercritical conditions, to remain in the non 
escaping regime for the flexible filament case. 


4- Non-escaping regime criteria. 

To avoid the presence of escaping filaments, one needs to have at equilibrium a negligible probability for filament 
size of the order of z*. This can best stated by comparing this probability to the (near) maximum value at i = z in 
supercritical conditions, giving 


n fi 

P{z) 


= a(z*,L)p[" ^ « 1 


where we have used Eq. (32 ). Taking the logarithm of both sides and using Eq. (371, one gets 

— p>W z * ( L ) + (z* - z) In pi <0 


(63) 


(64) 


The mean force potential W z * ( L ) is the reversible work to compress a grafted filament of size 2 * until it fits within 
the space limited by a hard wall at L. Using the approximate universal expression of the force in the weak bending 
limit Eq.(51l and treating the plateau value f bz * as constant over the whole compression interval, one gets 


-pfbz*C^ - !)£+ (| - 1 )^ ln ^i < °> 


Substituting Eq. (52) for /f, 2 * in Eq.(65l, the final expression of the non escaping regime condition reads 


'Ld 

pi < exp ( jy 


(65) 


( 66 ) 


which can be used either to limit pi at given L or to limit L at given pi. 

A comment about the evolution of size populations in the intermediate size window z < i < z* for any situation 


where condition (63) or equivalently (66l is met is in order. According to Eqs.(63), (641, one has 


P{i)/P{z) ~ exp ( (i — z) 


n 2 £ 

i a ip 

In Pi - 


4 d (i- l) 2 


Z < l 


(67) 


where we have again assumed the compressive force to be constant over the whole compression int erval (WpL) = 


( L c — L)f bi = [(* — 1 )d — L\TT 2 ( p /[4(i — l) 2 d]) and we have assumed L ~ (z — 1 )d according to Eq.(44l. The argument 


of the exponential is the product of a positive term (z — z) and the factor in square brackets where the first constant 
and positive term is dominated by the negative second term at the lowest i = z + 1 values as condition (|66| is met. 


According to Eq.(67l, P(i) must diverge as P(i) oc p\ when i grows to infinity. Therefore, the ratio in Eq.(67l must 
pass through a minimum (lower than unity) at some size i m in■ So if z* < i m i n (low pi value), the ratio P(i)/ P(z ) 
decreases monotonously over the relevant physical regime limited to z*, down to a small value required by Eq.(|66|). 
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FIG. 4. Single F-actin filament with l v = 5370d. L dependence of the ratio P(z*)/P(z) (in logarithmic scale) of the probability 
to have z* monomers over the probability to have z monomers. We report the ratio for two values of the free monomer 
reduced density pi = 1.7 (blue, continuous line) and pi = 2.5 (red, dashed line). The horizontal green line marks the value 
P(z*)/P(z) = 0.001 for which we consider that the bias induced by the constraint is negligibly small. The figure suggests 
maximal values of L to be L max /d ~ 89 for pi = 1.7 and Lmax/d ~ 70 for pi = 2.5 (see also Eq.( 


corresponding to P(z*)/P(z) = 1 level is reached for L/ = i p d/\npi which directly follows from Eq.( 661 . 


68 ). The red horizontal 


Otherwise, if i m in is located in the relevant z < i < z* regime (higher pi value), the criteria (63 ) implies that the 
ratio P(imin)/Pz must be even lower than P(z*)/P(z) so that kinetically, small filaments growing against the wall 
will see their size limited at values below i m i n . 

Eq. (66 ) predicts that the range of L values where the wall can effectively stop the bundle polymerization in 


supercritical conditions is limited by Li = \Jl p d/ In pi, namely Li ~ 100 and Li ~ 76 for pi = 1.7 and 2.5 respectively. 
In practice, for fixed pi > 1, we will limit the non escaping regime at the lower value L max which corresponds to 
P(z*)/P(z ) = 0.001, as illustrated in Figure [4] The dependence of L max upon pi, empirically established, is 


Lmax(pi) = Li(pi) - A(pi) = JIpd/ In Px - A(pi) 
A(pi) = 24.538 - 10.695 pi + 1.3965 p\ 


( 68 ) 

(69) 


where the A term has been fitted in the p\ range between 1.7 and 4. This relation provides L max — 89 and L max — 70 
for pi = 1.7 and 2.5 respectively, as seen in Fig|4j As a point of comparison, the measurement of the polymerisation 
force in an optical trap set up for what appears to be a single actin filament in supercritical conditions at pi = 1.7 
[5] (see Fig 4b of that reference), involves an elongation of L ss 200nm which corresponds to about L = 7Ad and is 
compatible with the present definition of the non-escaping regime. 

Condition (661 indicates that the concept of "non-escaping regime" is valid for flexible filaments only, since when 


l v —> oo the inequality is satisfied for any finite pi value. 


IV. THE STALLING FORCE AND ENERGY CONVERSION FOR F-ACTIN. 


A. The rigid living filament case. 


The living filament polymerization force F^(L, pi, T), shown for various cases in Figure [5l is derived in Eq. (601 
as the average of the L dependent force over an interval equal to the single monomer size d. Tn the limiting case of a 
rigid filament, any a(i,L) is a st ep function be ing unity as long as L > L C i (implying i < z(L)) and zero otherwise. 
Hence one has, using the Eq.(33l of D and Eq.(60) of D shl P, 


D 


shift! t\ _ 


(L) = D(L) = 


Pli^-Pl 1 ) 

1 - Pi 


(70) 
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FIG. 5. The polymerization force, averaged over an interval from L = nd to L = (n + 1 )d and reduced by ksT/d, is shown 
as a function of n for a single living filament at supercritical conditions pi = 1.7 (black curve) and pi = 2.5 (red curve). At 
each density, the flexible d-WLC case ( i p = 5370d, data points) is compared to the rigid limit case ( i v = oo, dashed lines) 
on the basis of Eq. ( |60| , which simplifies to Eq. (711 for the rigid case. The peculiar short L behavior, common to all curves, 
essentially reflects the boundary effect related to the imposition of a filament minimum size imin = 2. Beyond L > L (see text), 
the rigid filament average force goes to the asymptotic value In pi in agreement with Hill’s expression. Curves are deliberately 
interrupted in the Figure at L = Lm ax which is the upper range of the non escaping regime according to Eq.(68|. 


Using this relation, the average force in Eq. (|60| can be recast for the rigid filament case as 

Pi - 1 


dFfl 

k B T 


= In pi + In 


1 + 


Pi 


>(L) 


Pli 


(71) 


The correction (second) term in Eq.(711 is numerically important at small L only. For it to be of order e, one has to 
go beyond L given by 


L/d « I NT 


f ln(pi - 1) - ln(e) 




ln(pi) 


(72) 


which follows from Eq. (71) and from the link between z and L in Eq.(44). This boundary problem for rigid filaments 
is illustrated in Figure |5] where it can be o bse rved that the Hill’s result, lnpi, is indeed valid asymptotically beyond 
a value of L = L « 8d computed from Eq.(721 at pi = 2.5 for e = 0.001. 


Similarly in the rigid filament case, the size increment as the wall position L is displaced by d is given by combination 
ofEq.pland Eq.dTOi, 


Alfa(L,pi)) = 1 + pi—In 
op 1 


1 + 


Pi - 1 


= 1 + 


~z 7 

LPi - 1 


pu ( l) - 
(i-*) 
p\ z ~ 1 -1 


pi 


(73) 


Again, the correction to unity vanishes for large L/d (L > L where L is provided by Eq.(721). Figure [o] shows the 
increment becoming asymptotically unity for the rigid filament case as L increases. As commented by Hill[B] and 
shown in Figure [T] the ratio of the reversible work performed by the polymerization force to displace the wall by 
a distance d over the corresponding chemical free energy (pii — /x lc )A £fu used to increase the average length of the 
filament, goes to unity asymptotically for the rigid filament case. 
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FIG. 6. The increase Alfu(L, pi)) of the average length of a living filament pressing against a wall as the latter is moved from 
L = nd to L = (n + l)d under supercritical conditions specified by pi = 1.7 (black curve) or pi = 2.5 (red curve) is reported 
according to Eq.( |61| l. The rigid model case (t v = oo, dashed lines) (also given by Eq. m and the flexible case {l p = 5370d, 
data points) are compared. The behavior at short L (L < L) results from the imposition of a lower end boundary condition on 
filament length, namely i m in = 2. 



L/d 


FIG. 7. L dependence of the ratio of the reversible work of the polymerization force WL,L+d over the interval L £ [nd, (n+ l)d] 
and the corresponding chemical energy estimated as Alfu(L, pi) fee Tin pi. The results are shown for flexible filaments (£ p = 
5370 d) at two free monomers reduced densities, pi = 1.7 (black points) and pi = 2.5 (red points). For L > L, the observed 
central plateau value of unity reflects a perfect energy conversion. For the flexible filament case, the ratio starts decreasing 
progressively as L approaches the upper limit of the non-escaping regime. 
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B. The flexible living filament case adapted to F-actin. 


The effect of flexibility for a living grafted filament on the average polymerization force and on its size increment 
as the wall is displaced by the monom er s ize d needs to be investi gated f or supercritical conditions in the regime 

while the lower limit turns out to be in 


IIIB4 


L < L < L max . The higher limit, Eq.(68 >, was justified in section 
practice identical for flexible and rigid cases as illustrated in Figure |5J for two p\ values. Adopting e = 0.001, one has 
for pi = 1.7, L = 12d and L max = 89 d and for pj = 2.5, one has L = 7.9 and L max = 70. 

Considering the general expression for the average force Eq. (60), we first establish that in the relevant L regime 
and for the model of a WLC hitting a hard wall, the correction term is necessarily positive given the inequality 

D shl ft(L) > D{L) (flexible filaments) (74) 


This inequality basically follows from the property that a(i,L), as given by Eqs. (471, (501, is a monotonously 
decreasing function when L decreases, or equivalently when p increases (see Eq.(48l). This property is intuitively 
obvious and verified by visual inspection illustrative figures in ref. m■ To justify this on the basis of Eq.(|50|, we 


note that it is a sum of decaying exponentials in terms of fj but with alternating sign. Grouping terms in pairs, the 
even k — th term and the odd (fc+1) — th term, we obtain an absolutely converging series. The justification of Eq.(|74| 
follows from the property that, provided the sums of terms in D(L) and D sh ' l f t (L) converge sufficiently fast (before 
reaching i = z* , thus well within the non-escaping regime), the two expressions can be compared term by term. 
The strict inequality a(j + 1, L + d) > a(j, L ) for each pair of corre spo nding terms follows from the fact that, while 
L c — L is identical, the corresponding reduced compressions fj (Eq.(48l) is smaller for the £)«*”/*(£,) term, implying 
a larger value for a. In Figure [5j the average force Fff for a flexible actin living filament is shown as a function of 
L. The resulting curve in the regime L < L < L max is never very different from the rigid case, and thus from Hill’s 
prediction. The averaged force is slightly above the Hill’s plateau value by a marginal 1 — 2 % in the upper domain of 
the non escaping regime. This can be interpreted by noting that the a’s are related to the fraction of possible chain 
conformations for a chain of given number of monomer in presence of a rigid obstacle at distance L < L c . Intuitively 
this number should increase with the chain flexibility which is equivalent of taking longer chains L c + d at larger 
distance L + d for given £ p . 


The effect of flexibility on the filament size increment in the regime L > L gives, starting with Eq. (621, 


A lfu(L, pi) 


, „ d , 

1 + Piwt- In 
dpi 


D shi f\L,pi) 


D(L,p i) 

= i + [i s f h l ; ft (L,pi)-J fll (L 1 pi)) 


(75) 


where lfu(L,pi) is given by Eq.(34l and pi) is given by the same expression with probabilities P(i) = 


a(i,L)p\/D(Fpi) for i > z substituted by P(i) shl f t = a shl ^ t (i,L)p\/D shl ^ t (L, pi) (see also Eq.(60l). 

On Figured Alfu(L, pi) computed with Eq.(62) for the flexible case is again found to be close to the rigid limit. 
The value is however a few percents higher than unity in the upper part of the non escaping filaments domain, a 


logical result arising from the bending fluctuations of the filaments. Note that the approximate expression Eq. (75) 
(data not shown) gives identical results, except in the L < L domain. 

In Figure [7J the ratio of the reversible work of the polymerization force over the corresponding chemical free energy 
used to polymerize the living filament, goes also to unity for the flexible case at least in the central domain of the 
non escaping filament regime. At larger L, the ratio becomes lower than unity by a few percents, indicating that the 
conversion of chemical energy into work becomes affected by the flexible character of the filaments. Obviously, the 
situation quickly worsens if the filaments start to escape, which would happen with large probability if L gets larger 
than L max by 5 — 10 monomer units (Eqs.(68|69|). 


V. DISTRIBUTION OF FILAMENT SIZES PRESSING AGAINST A FIXED WALL FOR F-ACTIN. 


In this last section, we analyze the influence of flexibility on the equilibrium distribution of filament sizes when a 
living filament in supercritical conditions, is stopped by a normal hard wall. The size distribution, given by Eqs.(321, 
( [33| , takes in the rigid limit the form of a truncated growing exponential (P(i; L, pi) oc exp (i(ln pi)) for (i = 2, z(T/)) 
and P(i; L , pi) = 0 for (i > z(L )) and for (i = 1)). In order to avoid large fluctuations of the probabilities of hitting 
filaments as L varies over a monomer size distance d, we discuss results for the size distribution of flexible filaments in 
terms of an average over wall positions, as discussed earlier for the equilibrium polymerization force. Let Qk{L , pi) be 
the probability to have a filament of relative size k = i — z(L) with respect to the fixed wall position L. The average 
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< Qk >m computed as the average of Qk over the interval L/d € [n,n+ 1], is shown in Fig. [8] for several values of L 
covering the entire non-escaping regime. In Fig. [9] we show the average fraction ( Xo) n of filament sizes touching the 
wall, obtained as the cumulative sum of < Qk > n over the positive values of k. 



FIG. 8. Average normalized distribution < Qk >n of filament relative sizes k = i — z(L) at pi = 1.7, for wall position averaged 
over the interval [ nd , (n + l)d] shown for various n = L/d values within the non escaping regime. The exponential rise is 
observed for filaments sizes avoiding the wall, while the decay for filament sizes hitting the wall becomes increasingly sharper 
as L decreases. 


<dT 

c 


s 





FIG. 9. Fraction (xo) of filaments hitting the wall divided by lnpi plotted versus L for two supercritical densities pi = 1.7 
(black circles) and at pi = 2.5 (red squares). A clear trend (xo) oc InpiL 2 for all data in the non escaping regime is indicated 
by the green continuous line. 


These results show the most spectacular features of semi-flexible filaments with respect to rigid case. The (average) 
equilibrium polymerization force F ^ for filaments like F-actin in supercritical conditions and in the non escaping 
regime, is observed to remain essentially L independent and equal to the standard rigid filament result of Hill F^ = 
In pi. However, the way this force is produced by the living filament is highly L dependent and it is essentially 
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obtained as the product of two factors with inverse L dependencies 


T? av 


(xo) 


n 2 ksTip 

T L 2 


(76) 


as seen in Fig. [9] This first order expression means that the required force is produced at wall position L by recruiting 
a buckled filament of length L c ss L with a weight (xo) oc lnpiL 2 , while with weight (1 — {xo)), the filament does not 
contribute to the force on the wall, being it shorter than L. Indeed, filaments in contact with the wall are mostly in 
a compressed state corresponding to a force into the plateau region (see Fig. 2 ) therefore providing a force /& given 
by Eq. ( 52 ) in the force expression Eq. (51). This means for a wall located at L = 50 d that the plateau is reached 
as soon as ( L c — L) > O.lti which is most often the case (see Fig. [8}. In eq. (761 we also replaced 1 /L 2 ss 1/L 2 


Finally, it should be noticed that Eq. (761 also predicts that (xq) ~ i p 1 hence, in the rigid limit t p —>• oo, (xq) —t 0 


which demonstrates the impulsive character of the force when the rigid filament hits the wall during its Brownian 
fluctuations. 


VI. DISCUSSION AND CONCLUSIONS. 

The ability of actin filaments to sustain in supercritical conditions a compressive force has predominantly been 
justified with the aid of the rigid living filament model, effectively a one dimensional model, both for the single 
filament case and for bundles of parallel filaments [J1 IfiHTH]. The ID filaments which are fluctuating in length as a 
result of (de)polymerizing steps are hitting a fluctuating obstacle usually subject to load by producing instantaneous 
kicks which result into a time averaged force biasing the obstacle brownian motion. These brownian ratchet dynamical 
models do lead to an effective polymerizing force, compatible with Hill’s force expression at stalling, which satisfies 
a specific velocity-load relationship either for the single filament case [7] or for few bundle models differing mainly 
by the longitudinal disposition of filaments When Hill’s prediction for a bundle of Nf > 1 actin filaments 

has been found to fail in interpreting experimental data [5], the role of flexibility has been invoked by introducing an 
ad-hoc maximum value for the force that a bunble can exert and beyond which the bundle buckles. However even 
this ad-hoc extension of the ID model did not provide a satisfactory interpretation of the experiments. 

In the present paper and in its future extensions, we develop a Statistical Mechanics theory for a flexible filament 
model and we show that, quite generally, it leads to a systematic and continuous evolution of the filament behavior 
from a rigid rod character at short contour lengths to a pronounced flexible character at longer contour lengths, 
ending ultimately with the so called pushing catastrophe limit [171 118] . when the filament(s) has(ve) acquired a finite 
probability to grow unimpeded by the wall and escape laterally. While our approach will involve ultimately multi¬ 
filament bundles, moving obstacles under various loads thus implying non equilibrium situations, in this work we have 
focused on the already very rich phenomenology offered by the basic equilibrium properties of a single grafted filament 
in supercritical conditions, as it hits a fixed wall oriented normally to its grafting direction. We have incorporated 
the living character and the flexibility of the filament explicitly into a statistical mechanics approach based on the 
reactive grand canonical ensemble, dealing with the model of a discrete WLC hitting a hard wall. This formalism 
has been illustrated by the F-actin/free G-actin reacting mixture restricted to a single kind of actin monomer-ATP 
complexes. 

The results and the new phenomenology which emerges from the present work can be summarized as follows 

• We provide a statistical mechanics justification of the popular Hill’s expression for the single filament stalling 
force [B^. Rigorously, this expression corresponds to the average force exerted by the filament (defined as the 
ratio of the mechanical work over a finite displacement and the displacement itself), as the wall moves reversibly, 
under mechanical and chemical equilibrium, over a distance corresponding to one monomer size d. The Hill’s 
expression is found to be strictly valid only for the rigid filament case and we derive explicitly the correction 
terms for the semi-flexible case. The correction is positive and L-dependent for the model of a hard wall hit by 
a discrete WLC (the force is larger than for the rigid case) but these flexibility effects are only of the order of 
the percent when the experimental value of the actin persistence length is used. So we conclude that the L- 
independent Hill’s expression remains a very good approximation for the stalling force of semi-flexible filaments 
like actin. It should be stressed however that the force exerted by the filament on a fixed wall in the reactive 
grand canonical ensemble, that is the force which was integrated to get the work over a finite displacement of 
size d, shows large fluctuations around the mean. These fluctuations which decrease progressively in amplitude 
as L increases, find their origin in commensurability effects related to the degree of matching of the filament 
contour length, necessarily an integer number of monomer sizes d, with respect to the gap width L. These effects 
become less pronounced at large L as the amplitude of tip transverse fluctuations due to bending become more 
important. 
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• Like for the rigid case [TO], the equilibrium size distribution of the living flexible filament whose net polymerization 
is stopped in supercritical conditions by a wall at position L, starts as a growing exponential, as long as the 
filament size is too short for its set of fluctuating configurations to interact directly with the hard wall. Filament 
configurations larger than the slab gap have zero probability for the ID rigid model, while for flexible filaments 
the size distribution generally presents a fast decay which involves some finite but rapidly decreasing probability 
to get filament contour lengths larger than L. This rapid decay results from a filament bending work penalty 
which systematically exceeds the gain in chemical free energy, as a result of polymerization steps beyond the 
largest size z(L) of non touching filaments. Given the mentioned large oscillations in equilibrium properties as 
the gap width L is varied over sub monomer length scales, the filament size distribution properties are better 
discussed in terms of their d-averaged (average over a d window around the wall position L). The knowledge 
of these distributions (function of slab gap L) allows to adopt a quantitative definition for the limit of the 
non-escaping regime. In particular we require that the probability of a planar filament configuration of length 
7 tL/ 2, the minimum length to laterally escaping, be three orders of magnitude smaller than the probability 
of having filaments of length just below L. Our work provides the opportunity to establish more precisely the 
characteristics of the crossover towards the escaping regime, a point of high relevance in in-vitro experiments |U |5]. 


• At stalling, in the non escaping regime, the quasi L independence of the d-average force is produced by the 
fraction of filament configurations hitting the wall. The cumulative probability (xq) (L, p\) of the size distribution 
involving hitting filaments has been shown to increase like (xq) oc In pi L 2 / i p . This observation is compatible 
with the buckled filament state of the large majority of hitting filaments of the ensemble, each of them exerting 
adiabatically the classical plateau force expression f b = ^ksTj^ ss ^ksTj^ (here, adiabatic refers to the 
assumption that the life time of a given filament size is long with respect to the microscopic relaxation time of a 
fixed contour length filament). In this way the product (xo)fb is compatible with the L independent stalling force 
expression of Hill, which allows us to pinpoint a major distinction between rigid and flexible living filaments, 
a distinction established here for the case of a single filament at equilibrium but which will be relevant for 
multi-filament bundles at and outside equilibrium. For finite £ p , the L —> 0 limit of very short semi-flexible 
filaments leads to a contact probability (xq) —> 0 and a buckling force fb oo, just like in the case of rigid 
filaments (t p = oo) at arbitrary L . Hitting the obstacle takes the form of instantaneous kicks both for single 
filaments and multi-filament bundles. For flexible filaments of given £ p , the fraction of hitting filaments grows 
quadratically with L while the force of each (buckled) filament decreases quadratically with L. If we consider a 
dynamical trajectory of a single living flexible filament at equilibrium against a wall at distance L, the fraction 
of time the filament is in contact with the wall is finite together with the associated (buckling) exerted force. 
For bundles of Nf > 1 filaments at equilibrium, supposed to act independently, the force is produced by the 
permanent recruitment of a subset of (xo)Nf filaments pressing each with the buckling force fb(L), the subset 
of hitting filaments permuting continuously among the Nf equivalent filaments as the result of continuous 


(de)polymerization steps. Finally, as L approaches ~Jt p dj In p-\ under stalling conditions, which coincides with 
the upper limit of the non-escaping regime, the fraction (a,’o) should approach unity as the polymerizing force 
exerted by a filament cannot exceed fb- A more quantitative analysis of the limit is provided by eq. (691 which 
shows that the probability to get escaping filaments starts to be non negligible when (xq) approaches 0.5. 


A conjecture is possible when extending the criteria to observe the pushing catastrophe to the stationary situation 
of a wall moving at constant velocity, pushed by the polymerizing bundle of Nf > 1 filaments and subject to a load 
Fl = 'yF s tai smaller than the stalling value (7 < 1) [HI EH]. To keep a constant velocity of the wall, the bundle must 
exert a force equal to the load Fl (or sightly larger if solvent friction is considered), therefore the number of filaments 
N 0 = (x 0 )Nf needed to press on the wall should go as N 0 (L) ~ 7 F sta i/f b (L) ~ L 2 in order to compensate for the 
variation of the single filament force f b ~ 1/L 2 . If we assume that the non-escaping limit in stationary conditions 
(v > 0) would correspond to the recruitment of all Nf filaments (or a permanent contact with the wall for the single 
filament case Nf = 1), the maximum gap tolerated should be at least a factor 7 -1 / 2 larger than the limiting value 
for the non-escaping regime at stalling. 

It is illuminating to compare the above considerations with two reported experimental measurements of the actin 
polymerizing force. The experiments of Footer et al. [5] use an optical trap set up to measure the stalling force 
of a few actin filaments in supercritical conditions. In particular, they report in Figure 4b data corresponding to a 
polymerizing force of F sa 1 pN at reduced concentration p\ = 1.7 which, as they observe, corresponds to the stalling 
force of a single actin filament. The average filament length is « 180nm for the chosen optical trap. Eq.(76l applied 
to this case would imply a contact time fraction of (xg) ~ 0.25 with a force intensity of fb ~ 4 pN, the probed 
filament length being indeed lower than the limit L max = 240nm predicted by Eq. (69). This experiment was in fact 
dealing with a bundle of Nf = 8 filaments but surprisingly enough the stalling force of a single filament was effectively 
recorded. The issue here is still under debate but, according to our present work, to detect a force eight times larger 
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at the same reduced concentration in free monomers, a trap force constant 5 — 10 larger would be required in order 
to avoid laterally escaping filaments. 

The Demoulin et al. experiment [3] probes the force-velocity relationship for a set of actin bundles[3], implying a 
total of Nf ss 130 filaments at pi ss 3, pressing together against a bead. While the stalling force in this case is around 
200piV, the bead is subject to load forces ranging from a few pN (largest velocity probed) up to lOOpiV covering a 
range 0.02 < 7 < 0.5 for the load over stalling forces ratio. Looking at figure 2 in ref.[3], if we take a typical length of 
200 nm for the actin filaments beyond their lateral connection by fascin bridges, the number of filaments at contact 
able to press on the obstacle bead with a buckling force of « 4 pN should lie between 1 and 25 over the explored force 
range. Further, considering the results for the longest filaments (~ 400nm) at Fl = S.dpN, our criterium above for 
stationary non-escaping conditions is still justified since L max (pi) = 65 d = Ylbnm at p\ = 3, and y^ 1 / 2 k, 7 in these 
conditions leading to a maximum length of the stationary non-escaping regime of L max / 7 1 / 2 = 1225nm, still larger 
than the probed bundle length. 

Further consequences of filament semi-flexible character on actin bundles at and outside equilibrium will be discussed 
in future publications ESI US]. 
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Appendix A: Some considerations on (de)polymerization rates. 


Phenomenological kinetic rate constants kon ^ and k l of j: are usually associated to the reactions described by Eqs.fij). 
In terms of such kinetic constants the equilibrium micro-reversibility conditions, i.e. the equality between the number 
of polymerizations of filaments of size (* — 1 ) to size i and the number of de-polymerizations of filaments of size i to 
size (* — 1 ) per unit of time, are written as 




which implies, according to Eq.(41) and thermodynamics, the link 

Ki = kt 1] /K ff 


Eq.( Al I is often written equivalently as 


U [i _ l) P{i-l) = W i P{i) 


(Al) 


(A2) 


(A3) 


where U^_ 1 ) = kon ^ Pi and W t = are (de-)polymerisation rates, respectively. Using again Eq.(|4l|, their ratio is 


U (j- 1 ) 

W, 


P(i) 




P(i-l) Pl a(i-1,L) 


(A4) 


The (de)polymerization rates for filament ends in bulk, denoted by Uq and Wq with Uq/Wq = pi, are valid for our 
grafted filaments, as long as they do not interact with the wall. When pi = 1 (that is the free monomer density is 
critical pi = pi c = \/Kq), one has Uq = Wq so that the rate of polymerization and rate of depolymerization are equal 
and the distribution should be uniform in the short filaments region (a = 1). We are interested to the supercritical 
regime pi > 1 and Uq > Wq for which the distribution in the same short filament region is a growing exponential. 
For (de)polymerization reactions implying filaments hitting the wall, the rates satisfying Eq.(A4l are often chosen in 
applications assuming that the rates of depolymerisation are not affected by the presence of the wall 0 nans 123, 
namely 


Wi = W 0 

a(i,L) 

u i - 1 — —r- —TT 3 
a(i — 1 , L) 


U 0 


(A5) 

(A 6 ) 
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